Exact Duality Relations in Correlated Electron Systems 
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Using gauge transformations on electron bond operators, we derive exact duality relations between 
various order parameters for correlated electron systems. Applying these transformations, we find 
two duality relations in the generalized two-leg Hubbard ladder at arbitrary filling. The relations 
show that unconventional density-wave orders such as staggered flux or circulating spin current 
are dual to conventional density-wave orders and there are direct mappings between dual phases. 
Several exact results on the phase diagram are also concluded. 
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Unconventional density-wave orders such as stag- 
gered flux [which is equivalently called as d-density 
wave (DDW) or orbital antiferromagnets] and circu- 
lating spin current were first proposed in the context 
of excitonic insulatorsUj and later discussed in high- Tc 
superconductors|3,IEHj but the appearance of these or- 
ders was not established at that time. Recently, several 
experimental results have led to a resurgence of interest 
in the possibility of the exotic orders. A circulating spin- 
current state'5'1 and DDW statejfij were proposed as an 
origin of hidden order in the heavy-fermion compound 
URu2Si2. A DDW state has also been discussed to ap- 
pear in underdoped region of high- Tc superconductors 'Ti], 
where a pseudogap was observed 8] , and in the quasi- two- 
dimensional organic conductor Q!-(ET)2KHg(SCN)4|^. 

From a theoretical viewpoint, the possibility of the un- 
conventionally ordered phases in microscopic models of 
correlated electrons has been the focus of interest. In 
particular, the two-leg ladder system is attracting atten- 
tion as a minimal model for showing the exotic orders. 
Though the two-leg Hubbard and t-J ladders do not show 
any order [lol Ull . both strong coupling and weak cou- 
pling analyses revealed that the generalized two-leg Hub- 
bard ladder including nearest-neighbor interactions ex- 
hibits various ordered phases at half-filling fislll^lislll^ 
and doped cases 0,0]. Large-scale numerical calcula- 
tions also reported the DDW phase [lH at half- filling and 
charge-density- wave (CDW)'l^ and DDW[l3 phases for 
less than half- filling. 

In this Letter, introducing gauge transformations on 
electron bond operators, we derive two exact duality re- 
lations between various exotic phases in correlated elec- 
tron systems. One duality relates conventional density 
waves to unconventional density waves (or currents) and 
the other relates charge-density degrees of freedom to 
spin-current ones. These transformations also show du- 
ality relations between s- and d-wave superconductivity 
(SSC and DSC), as well as between s- and d-wave Mott- 
insulating (S-Mott and D-Mott) states. Applying these 
transformations, we show two exact duality relations in 
the generalized two-leg Hubbard ladder (or similarly the 



Hubbard model with doubly degenerate orbitals). The 
duality relations among various phases are summarized in 
Fig.^ where we use density- wave terminology^^ . Start- 
ing from previously established phases, we can reveal the 
appearance of various exotic phases with and without 
spin symmetry breaking. The transformations give pa- 
rameter mapping between dual phases and fix boundaries 
between them on self-dual lines. These exact relations en- 
able us to clarify the nature of the phase diagram and find 
the parameter space of unconventionally ordered phases. 

Before discussing the duality transformations, we list 
operators and order parameters considered in this work. 
We use ladder geometry in the following for simplicity, 
but one can extend to other lattices adopting appropriate 
bonds. Regarding charge degrees of freedom, we consider 
CDW, staggered-flux [or equivalently DDW], staggered 
dimer [or p-density-wave (PDW)], and diagonal current 
[or /-density-wave (FDW)] operators on rungs (or pla- 
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FIG. 1: Duality relations (I) and (II) between various states; 
(a) CDW, DDW, SDW, and DSDW states, (b) PDW, FDW, 
PSDW, and FSDW states, and (c) s- and d-wave Mott- 
insulating and superconducting states. 
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quettes) given by 
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CFDw(i) = ^5I[(4 + 1,1:<tCj.2,o- -c] + 1,2,<tCj,1,^) 

-H.c], 



(4) 



respectively, where c]^^ ,^ (cj>,ct) (m = 1,2 and cr =t,i) 
denotes an electron creation (annihilation) operator on 

the /ith site of the jth rung and nj,p.(T = c^j f_i a'^j,t^,<^- 
We note that the CDW and PDW orders are kinds of 
density waves while the DDW and FDW orders finite lo- 
cal currents with broken time-reversal symmetry. In the 
similar way, we consider operators in the spin sector. In- 
serting the Pauli matrix cr^ ^ into the right-hand sides 
of Eqs. Q-Q, one defines spin-density wave (SDW) op- 
erator OsDW [= h^a^^a.ai^jA.cr - nj,2,a)], circulating 
spin current [or d-spin-density wave (DSDW)] operator 

Cdsdw [= 5 J2a '^l,a{^].i,<y^oa.'y " H-C-)], p-spin-density- 
wave (PSDW) operator Cpsow, and /-spin-density wave 
(FSDW) operator Ofsdw as the counterparts of the 
CDW, DDW, PDW, and FDW operators, respectively. 
Order parameters on the ladder are given by 

OA(<7) = i-'^OA(j)exp(zgj), 
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where L is the number of rungs. Furthermore, we con- 
sider the d- and s-wave pairing operators, 

CdscO') = ;^(ci,i,TCj,2a - CjM,iCi,2,T), (5) 

CsscO') = ;^(ci4,TCj,i,x +Cj,2,TCj,2a), (6) 

which characterize the DSC and SSC states. We also 
consider the D- and S-Mott states at half-filling, 

|D-Mott) - n^DSc(j)|0>, |S-Mott) = n^^ssc(j)lO), 



where |0) is the vacuum of the electron operator Cj,^,ct- 

Let us move to the introduction of duality transfor- 
mations on electron operators. They are given by gauge 
transformations on the bonding and antibonding opera- 
tors dj^i^o- = (cj,i,cr±Cj_2,dr)/'\/2. We consider two duality 
transformations and their combination in the following. 

Duality I: density and current. — Let us consider a 
gauge transformation of antibonding operators given by 



the unitary operator Ui{6) = Y\j ^Qy^vi—iOd'^ 
When 9 = 7r/2, this operator yields the duality transfor- 
mation (I) as 



(7) 

for CT =t, J,. In terms of the electron operators Cj,^,a, this 
transformation is written as Cj^i^a = ^^Cj,i^a ^^Cj,2,a 
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One finds that the 



and Cj^2,a — ~2^Cj^i^cr 

duality transformation relates the density- wave operators 
to the current ones as follows: 



(8) 



Thus, current phases, i.e., DDW, DSDW, FDW, and 
FSDW phases, turn out to be dual to density- wave 
phases, i.e., CDW, SDW, PDW, and PSDW phases, re- 
spectively. In other words, this transformation changes 
the bond parity and hence s- (p-) wave orders are con- 
verted to d- (/-) wave ones. The pairing operators and 
the Mott states are transformed in the similar way, 
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We note that the unitary operator Ui{9) gives a con- 
tinuous transformation between dual operators, e.g., 
Ui{0)OcBwUiie)-^ = OcDwcosS - OoDwsine. 

Duality 11: density and spin current. — Next, we con- 
sider a gauge transformation given by the unitary op- 
erator Uii{9) = njexphj6'(d]_+jdj,+,T + d''-_^^dj-^i)]. 
This operator with 9 — t:/2 gives the duality transfor- 
mation (II) as 



d_j,+,l - idj^+^^, 
dj-A = '^j -,T> 



(9) 



This transformation converts the operators as 



CcDW — Cdsdw, 
Cddw — — OsDW, 
CpDW = C'fsdw, 

OfDW = — C'PSDW, 



Cdsdw = ^Ccdw, 
CsDW = Cddw, 
CpsDw = — CpDw, (10) 
CpsDw — Ofdw- 



One can see that density waves of charges are trans- 
formed into spin currents while density waves of spins 
into charge currents. This transformation thus exchanges 
density and current as well as spin and charge degrees 
of freedom. Similarly to Ui{9), the unitary operator 
Uii{9) gives a continuous transformation between dual 
operators such as Uii{9)OcdwUii{9)~^ = Ocdw cos 61 -I- 
Cdsdw sin 9. We note that under this transformation the 
pairing operators and Mott states are invariant except for 
phase factors. 
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Duality between spin and charge. — Combining the du- 
ality transformations |7J and © , we obtain the transfor- 
mation 



2i: + ,T 



— dj 



= —id 
— idi 



(11) 



From the duaUty relations and Hl()|l . it follows that 
this transformation exchanges spin and charge degrees of 
freedom, e.g., density- wave operators are transformed to 
spin-density-wave operators. 

To elucidate efficiency of these transformations, we ap- 
ply them to the generalized two-leg Hubbard ladder with 
intra-rung couplings. The Hamiltonian is defined by 



H 



+ipair(c],l,tc] ;^^^Cj,24Cj^2,T + H.C.) , 



(12) 



where n 



T,a %,A','T and Sj, , 



^J2a,cr' ^j.ti.a'^'^.<^"^J,t^,cr' ■ When t±_ ^ 0, this Hamilto- 
nian includes the Hubbard model with doubly degenerate 
orbitalS|l^. We first apply the transformation (I) [Eq. 
^] to the Hamiltonian. It can be shown that this 
transformation maps the model onto the same Hubbard 
model with different coupling parameters. The total 
charge density, the total magnetization, and the kinetic 
energy terms are invariant under the transformation. 
However, the coupling terms are mixed by the trans- 
formation. Here, it is useful to rewrite the couplings in 
terms of bonding and antibonding operators as 



(d) (d) 



(13) 



Q = ± 



a=± 



where 4\2,<t = d]^„ „d.j^a,a, and 

A = ([/ - + ipair)/2 - 
{U-V^- ipair)/2 - 

Vi_ + Ji/4, 
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r)/2 



{U + V^~ Vir)/2 + (2 J± 



Jl)/8. 



It is easily found in Eq. that the transformation 

O changes only the sign of the A-term, but keeps the 
rest of terms invariant. This leads to an exact duality 



relation in the two-leg Hubbard ladder. The model with 
a parameter A is dual to the model with —A and the 
model is self-dual in the space A — 0, i.e.. 



U-V±+ tpair + (2 J± + Jl)/4 - 0. 



(14) 



The mapping of the original coupling parameters is as 
follows: 



u = i(;7 + v^i-Vir)-^(2J± + Jl), 

Vl_ = i(;7 + 3Vl+ipair) + ^(2J± + JI), 



J± = -U + V^-tp,i, + -{2J^- Ji), (15) 

Ji = -C/+l/±-tpair-^(2Ji-3Ji), 
ipair - + V±+tp,i,)- 1(2.1^ + J^). 

Note that spin isotropy {J± = J^) is conserved through 
this parameter mapping as J± = J^. From the map- 
ping (|15|l and the duality relation ((SJ, one can conclude 
that if a density-wave order, e.g., CDW or PDW order, 
appears in a certain parameter region, a dual current or- 
der, i.e., DDW or FDW order, respectively, exists in a 
corresponding dual region of the parameter space. Be- 
cause of this duality relation, all phase boundaries must 
be symmetric with respect to the self-dual space. Indeed, 
the transition line between the CDW and DDW phases 
at half-fil ling derived in the weak- and strong-coupling 
limits m M, Ei 

coincides with the self-dual line p4|l . 
We stress that our result holds in general cases, regard- 
less of the coupling strength, filling, and system size. 

Next we discuss the duality transformation (II) [Eq. 
©], which gives another parameter mapping. Similarly 
to the transformation (I), the total charge density, the 
total magnetization, and the kinetic energy terms are 
invariant under the transformation but the coupling 
terms are mixed up. From Eq. 113|) . one can see that the 
transformation changes only the sign of the i3-term. 
We thus find that the present model has another duality: 
The model with a parameter B is dual to the model with 
—B under the transformation ^ and the model is self- 
dual in the space 5 = 0, i.e., 

U-V^- ipair - (2Ji - Ji)/4 = 0. (16) 

The transformed coupling parameters are given by 

2* 



U = -{u + V^+tp,i,) + -i2J^-Ji), 



V± = i([/ + 3l^i-ipair)-^(2J±- Ji), 



J± = C/-y±-ipair + 7(2J± + JI), (17) 



JI - -C/+l^i+tpair + |(2J±+3JI), 
ipair = ^{U~V±+tp,i,)-^{2J±-Jl) 
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This duality relation leads to the conclusion that, if CDW 
or DDW order, for example, exists in a certain parameter 
region, spin current (DSDW) or SDW order exists in a 
dual region, respectively. Note that even if we start from 
a spin isotropic model the dual model is spin anisotropic 
and hence spin symmetry breaking can occur in the dual 
model. The phase diagram must be symmetric with 
the self-dual space (fTH|l and the transitions between dual 
phases, if ever, locate exactly on the self-dual space. 

The spin-charge duality relation associated with the 
transformation can be obtained from the combina- 
tion of two transformations discussed above. Thereby 
one can arrive at a direct duality relation between orders 
in the charge and spin sectors. The self-dual space is the 
intersection of Eqs. and JTHl), i.e., U-V± + \JI = 
and 5 J± -I- tpair = 0. 

In the self-dual spaces H14|) and H16|l . the model Hamil- 
tonian is invariant under the continuous transformations 
with the unitary operators Ui{9) and Uii{9), respec- 
tively, and hence U(l) symmetric. [This can be easily 
seen in Eq. itHni]. Because of the U(l) symmetry, a rig- 
orous theorem 20!| concludes that the dual orders that 
are continuously transformed to each other disa ppe ar on 
the self-dual models in one dimension, in generalj2lj, and 
the phase transition between the dual phases, if ever, is 
of second order. The A- and _B-terms in the interaction 
(|13|l serve as symmetry-breaking perturbations for the 
U(l) symmetry of Ui{d) and Uii{6), respectively. When 
a symmetry-breaking perturbation is relevant, it induces 
order and a gap determining the criticality of the phase 
transition. 

In summary, we have developed duality transforma- 
tions for electron systems, which lead to duality rela- 
tions between various exotic phases shown in Fig. ^ 
These duality relations have an analogy with the spin- 
chirality duality transformation we and co-workers intro- 
duced for the spin ladder 'i^, E^, but are applicable to 
a much wider variety of systems. The transformations 
clarify for the generalized two-leg Hubbard ladder that 
the stability of unconventional density-wave orders such 
as staggered flux (DDW) and circulating spin current 
(DSDW) is equal to that of conventional density-wave 
orders in the dual spaces. Recent large-scale numerical 
analyses ^2 Il8«l reported the appearance of CDW and 
DDW phases under doping. Application of the dual- 
ity relations to these results immediately concludes that 
SDW and DSDW phases stably exist under doping in 
dual parameter spaces. Further arguments and possi- 
ble applications of the duality transformations will be 
reported elsewhere. 

Finally we stress that the duality transformations (0 
and JHJ can be applied to a variety of systems in which 
bond degrees of freedom are important, since they are 
based only on simple gauge transformations of bond op- 
erators. Applications to double- layer systems, and two- 



and three-dimensional systems with orbital degeneracy 
are of interest. 
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